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ABSTRACT 
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I.    INTRODOCTION 

A-       EAfiTHQOAKE   PREDICTION    THEOfiT 

The  desirability  of  being  able  to  predict  when  and  where 
an  earthquake  will  occur  becomes  immediately  apparent  in 
light  of  the  recent  devastation  that  occurred  in  Mexico 
City.  The  potential  for  death  if  a  catastrophic  earthquake 
should  cccur  along  the  San  Andreas  fault  is  estimated  to  be 
between      3,000    and      13,000    lost     lives.  The    lower      figure 

applies  for  a  quake  that  hits  at  0230  and  the  higher  figure 
applies  for  a  quake  that  occurs  at  1630  on  a  weekday.  These 
predictions  are  based  on  the  facts  that  wooden,  single 
family  homes  are  much  more  resistant  to  earthquaKe  struc- 
tural damage  than  are  many  of  the  older,  unreinforced 
masonry  buildings.  In  May,  1983,  30  of  the  40  unreinforced 
masonry  structures  in  Coalinga,  Ca-,  either  partially  or 
totally  collapsed.  In  Los  Angeles,  8,000  similar  buildings 
exist.  The  probability  of  a  magnitude  8.3  earthquake  occur- 
ring along  the  southern  San  Andreas  fault  within  the  next  30 
years   is   estimated   to   be    50    percent-      [ Ref .    1:p-    35] 

The  primary  cause  of  earthquakes  is  the  movement  of  the 
litaosphere,  the  outer  layer  of  the  Earth  consisting  of 
twelve  riqid  plates  riding  on  the  supporting  mantle.  This 
movement  rearranges  continents  and  forms  and  destroys  botn 
mountains   and    oceans-  The    driving    energy    for      this    poten- 

tially large  scale  destructive  mechanism  is  the  convective 
circulation   of      the    mantle.  The    mantle      consists    of      hot, 

solid  rock-  In  time  measured  on  a  geological  scale,  the 
mantle  is  sufficiently  hot  that  it  can  flow  and  be  easily 
deformed-  Hot  magma  from  the  mantle  walls  up  between  sepa- 
rating plates  at  midocean  ridges  forming  new  lithospnere. 
The  movement  of  the  cold,  brittle  plates  is  primarily  along 
boundary   faults   and    is  controlled    by    friction.       Earthquakes 
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occur      when      these      plates      suddenly      slip      and      elastically 
rebound.      £Ref.    2] 

Figures  1.1  and  1.2  [Eef.  l:p.  38]  show  the  typical 
cycles  of  earthquakes  as  explained  by  the  theory  of  plate 
tectonics.  In  figure  1.1  the  hot  mantle  10-15  km  deep  is 
plastically  deforming  and  stress  is  accumulating  at  asperi- 
ties, or  lock  points,  in  the  cooler  crust.  At  the  moment  of 
rupture  an  asperity  breaks  and  crustal  rupture  spreads 
causing  earthquakes-  The  slippage  is  halted  when  it  encoun- 
ters another  asperity  in  approximately  100  seconds. 
Aftershocks  may  occur  for  the  next  several  months. 
Additional  slippage  at  or  near  the  surface  may  or  may  not 
occur.  In  Figure  1. 2  a  second  major  earthquake  has  occurred 
50  or  100  years  after  the  first.  The  cycle  begins  again 
after  about  200  years  when  the  deep  crust  displacement 
imposes   new   stress   on   the   upper   brittle   crust. 

To  date,  earthquake  prediction  has  been  on  a  long-term 
basis  normally  expressed  as  a  probability  of  quake  occur- 
rence during  the  coming  scores  of  years.  What  has  been 
lacking  is  a  reliable,  short-term  prediction  process. 
Central  to  any  short-term  technique  is  geologic  phenomena 
during  the  period  immediately  preceding  the  main  earthquake. 
Before  the  arrival  of  a  major  shock,  swarms  of  many  small 
earthquakes  are  often  noted  or  sometimes  accelerated  defor- 
mation of  the  crust  is  observed.  However,  neither  one  of 
these  two  predictors  is  reliably  present  in  all  cases. 
[Ref.    1:p.    38] 

Systematically  collected  data  on  the  generation  of  elec- 
tromagnetic energy  prior  to,  during,  or  after  an  earthquake 
currently  does  not  exist.  Evidence  of  radio  frequency  (RF) 
emissions  as  a  precursor  to  quake  activity  appears  impul- 
sively suggesting  these  waves  are  generated  over  a  broad 
spectrum  of  frequencies  [Ref.  3].  Citizen  band  and  amateur 
radio  operators  in  the  Hollister,  Ca. ,  area  reported  an 
"increase      in    background      noise      level    preceding      earthquake 


1  1 


^.^1. 


"^i'-V,,  ,    ^  ^  ''H 


,.-:::iJ 


Figure  1.1   First  Portion  of  an  Earthquake  Cycle 
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Figure  1.2    Second  Portion  of  an  Earthquake  Cycle 

activity"  at  27  MHz  and  at  the  six  and  two  meter  bands. 
This  increase  in  noise  was  on  the  order  of  10  dB  and 
preceded  a  quake  by  12  to  24  hours-  At  one  hour  to  1 5 
minutes  before  a  quake,  noise  levels  returned  to  normal. 
The  appearance  of  earthqaake  lights  [Ref.  4]  and  reception 
of  18  MHz  radio  noise  [Ref.  5]  are  further  examples  of  a 
possible   connection   between   RF   energy   and   earthquakes- 
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Various  theories  explaining  the  actual  mechanism  of  SF 
generation  at  the  fault  boundary  include  impact  of 
dissimilar  metals,  shock  induced  polarization  and  fracture 
of  quartz-bearing  rock  [Ref-  6].  Because  of  the  strong 
implication  that  earthquakes  generate  electromagnetic 
energy,  this  thesis  will  present  a  satellite  design  for  a 
space-based  system  as  a  means  of  efficient  quake  detection 
and    rapid  data    dissemination. 

B.       SFACE-BASED    EARTHQOAKE    DETECTION 

The  exploitation  of  outer  space  as  a  tool  in  earthquake 
prediction  based  on  electromagnetic  detection  can  be  divided 
into  two  major  scenarios.  The  first  is  installing  ground- 
based  sensors  that  detect  the  quake  generated  fiF  energy  and 
then    relay      data    to   an      orbiting   satellite.  The    satellite 

would  function  as  a  central  collection  point  and  could 
incorporate  features  such  as  data  processing  and  earthquake 
pattern  recognition  with  automatic  alarm  transmission  to  key 
Earth  stations.  The  number  of  ground-based  sensors  to 
provide  adequate  coverage  can  not  be  reliably  estimated 
because  the  strength  of  the  RF  energy  is  not  known  but  is 
being  researched  as  a  portion  of  this  thesis.  As  the 
strength  increases,  the  distance  between  sensors  can 
increase   and    their   number    decrease. 

The  second  scenario  is  direct  detection  of  earthquake 
generated  £F  energy  by  an  orbiting  satellite.  Ail  the 
previously    mentioned      features    could      be    included.  Direct 

satellite  detection  would  require  a  sufficiently  strong  RF 
signal  to  overcome  the  power  spreading  loss  incurred  from 
the  large  distance  traveled.  The  major  advantages  over  the 
first  scenario  are  worldwide  coverage  witn  elimination  of 
the  need  for  ground- based  RF  sensors.  On  an  international 
level,  frequencies  could  be  allocated  to  quake  detection  to 
minimize  possible  man-made  interference  with  the  RF  signal. 
Because  of  tne  advantages,  this  second  scenario  will  be  the 
basis    for   the    proposed   satellite   design. 
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C.       FCCUS    OF    STUDY 

A   proposed      design  for    a      space   based   sensor      to    predict 
and    detect      earthquakes   is    presented-  A    free      standing   EF 

noise  measurement  and  recording  system  is  described  to 
research  the  correlation  between  earthquakes  and  increased 
background  electromagnetic  noise  at  30.45  MHz  and  150.75 
MHz.  This  thesis  continues  the  research,  construction  and 
installation   of  equipment   described   in   Reference   6. 
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II.     BACKGBOOND 

A.       OfiBIIAL    PAEAMETEfiS 

1 .      Space   Coordinate   System 

In  order  to  understand  the  bacieground  information 
that  will  be  used  in  subsequent  portions  of  this  thesis,  it 
is  necessary  to  first  establish  the  coordinate  system  that 
describes  a  satellite's  position  in  space.  An  ephemeris  is 
a  tabulation  of  the  predicted  or  observed  positions  occupied 
by    the   satellite. 

The  coordinate  system  to  be  used  is  shown  in  Figure 
2.1    [  Bef -    7]-  The    large,         flat    rectangular      plane  is      an 

extension  of  the  plane  of  the  equator  and  is  called  the 
equatorial  plane.  Superimposed  on  the  equatorial  plane  is 
the  Cartesian  coordinate  system  indicated  by  I,  J  and  K  with 
the  origin  at  the  Earth's  center-  These  coordinates  do  not 
rotate  but  are  fixed  in  space  with  the  I  axis  pointing  in 
the    direction    of    the    vernal    equinox. 

The  satellite's  motion  about  the  Earth  is  on  the 
circumference  of  the  orbital  plane.  An  elliptical  orbit  is 
shown  but  a  satellite  can  also  have  a  circular  orbit  around 
the  Earth.  Tns  angular  momentum  or  spin  vector  (h)  for  tae 
satellite  orbit  is  determined  by  the  right-hand  rule  for  a 
south  to  north  equator  crossing  and  is  perpendicular  to  the 
orbital  plane.  Inclination  (i)  is  measured  in  degrees  from 
K  to  the  spin  vector  h.  By  definition,  inclination  can  vary 
from    a    minimum   of    zero   degrees    to   a    maximum    of    180    degrees. 

The   eccentricity    vector    (e)  starts    from   the   origin 

and  points  in  the  direction  of  perigee  or  the  periapsis 
dicection.  Perigee   is      the    satellite's      closest    point      of 

approach  to  the  Earth.  The  line  of  nodes  ^nj  is  the 
intersection  of  the  equatorial  plane  and  the  satellite's 
orbital  plane.  The  angle  between  the  line  of  nodes  and  the 
eccentricity    vector    is      called    the    arguement   of      perigee    and 

16 


Figure  2.1    Classical  Orbital  Elements 
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is  measured  in  the  direction  of  the  satellite  motion. 
Right  ascension  of  the  line  of  nodes  (n)  is  the  angle 
measured  in  the  equatorial  plane  from  the  vernal  equinox 
direction  to  the  line  of  nodes  "n.  The  epoch  angle  {..^  )  is 
the  angle  measured  from  e"  to  the  vector  from  the  center  of 
the  Earth  to  the  satellite  (r„  ).  The  range  of  the  epoch 
angle  is  from  zero  degrees  to  360  degrees. 
2.      Equations   for  an    Ellipse 

Many  satellites  orbit  the  Earth  in  an  elliptical 
manner  or  use  portions  of  an  elliptical  orbit  daring  an 
orbital  transfer  from  one  circular  orbit  to  another  circular 
orbit  of  a  different  radius.  Because  the  equations  of  an 
ellipse  are  an  integral  portion  of  satellite  orbitology,  the 
pertinent  relationships  will  be  presented.  Figure  2.2  shows 
the  classical  form  of  an  ellipse  with  the  small  circle 
representing    the    Earth  at   one    focus. 

The  semi-major  axis  (a)  is  the  horizontal  length 
measured  from  the  center  of  the  ellipse  to  the  outer  edge. 
The  semi-minor  axis  (b)  is  the  corresponding  vertical 
length.  The  eccentricity  (e)  is  a  measure  of  departure  from 
a   cirlce   and   is  given  by   Equation   2.1 

e    =   SQIiT(1-{b/a)2)  (2.  1) 

Eccentricity  can  vary  from  a  value  of  zero  indicating  b  =  a 
for  a  circle  to  a  value  of  slightly  less  than  one  if  b 
approaches    zero. 

The  radius  of  apogee  (Rq)  is  the  distance  from  the 
center  of  the  Earth  to  the  point  of  farthest  distance  from 
the  Earth.  The  radius  of  perigee  (Rp)  is  measured  from  the 
center  of  the  Earth  to  the  satellite's  closest  point  of 
approach  to  the  Earth.  The  relationship  between  Rq  and  e 
can    be   expressed    as 

Rq    =   a  (1+e)  (2.2) 
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Figure    2.2        An   Ellipse  Around    the   Earth 
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A   similar   formula    holds  for    Rp   and   e, 


Rp    =   a{1-e)  (2.3) 


B.       ORBIT    EQOATIOHS 

1-      Circular    and    Elliptical   Orbit   Velocities 

Two  celestial  bodies,  one  in  orbit  around  the  other, 
comprise  a  two-body  system  and  must  obey  Newton's  law  of 
universal   gravitation   given    by 

F    =   GMm  (2.4) 

r2 

Referring  to  Figure  2.3,  F  is  the  force  of  attraction  due  to 
gravity  between  masses  M  and  m  separated  by  a  distance  r- 
The  constant  G  is  the  universal  gravitational  constant  and 
its  value  is  6.6  7  X  10-ii  m^/kq  sec^.  For  any  body  of  aass 
ra  in  a  circular  orbit  with  velocity  v  and  at  a  distance  r, 
the    centripital    force    is   given    by 

F    =   mv2  (2-5) 

r 

The  centripital  force  must  equal  the  gravitational  force  or 
the    orbiting   body    will  not    be    in   an    equilibrium      orbit. 

G3m      =      mv2  (2.6) 

r2  r 

Cancelling  the  common  terms  of  m,  one  of  the  r's  and  solving 
for    v    yields 

V    =    SQfiT(GM/r)  (2.7) 

This  V  is  the  linear  velocity  of  mass  m  in  a  circular  orbit 
about  mass  M  at  a  radial  distance  r.    The  time  required  for 
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m  to  complete  one  circular  revolution  (Tj.)  is  found  by 
dividing  the  distance  traveled  (2r7r)  by  the  velocity.  This 
gives 

Tg     =   2r7r  (2.3) 

V 

=   2r7r  (2.9) 

SQRT(GM/r) 

=   2  TT  r3/2  (2.  10) 

SQET(GM) 

Kepler's  third  law  of  planetary  motion  relates  the 
motion  of  one  planet  to  another.  This  relationship  is  inde- 
pendent of  the  shape  of  their  orbits  and  requires  that  the 
sguares  of  the  periods  or  the  time  to  complete  one  revolu- 
tion vary  as  the  cubes  of  the  semi-major  axes.  Then  for  a 
in  an  elliptical  orbit  about  M  with  a  semi-major  azis  of  a, 
the    time   for    one    revolution     (Tg)    is    given    by 

Te    =   27ra3/2  (2.  11) 

3QRT(GM) 

Tha    elliptical    orbit    velocity     (Vg )    is    given    by 

Vg     =    SQRT[Gi1  (2/r-1/a)  ]  (2.12) 

2-      Orbital    Transfers 

It  is  often  desired  to  change  the  orait  of  a  satel- 
lite to  place  it  either  in  a  more  advantageous  relative 
position  to  other  satellites  or  to  alter  its  relative  posi- 
tion to  the  Eartn  and  allow  it  better  access  to  Earth  gener- 
ated data.  To  perform  an  orbit  modification  the  satellite 
must  expend  energy,  normally  in  the  form  of  burning  fuel 
from    an    on-board    rocket   or      releasing    compressed    gas    tnrougn 
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Figure  2-3    Gravitational  Force  of  Attraction 

thrusters.  Rockets  are  used  when  a  major  change  is  desired 
because  a  large  force  can  be  quickly  developed  while 
compressed  hydrazine  gas  powered  thrusters  are  normally  used 
for  station  keeping  or  very  minor  orbital  adjustments. 

The  two  major  factors  that  influence  exactly  how  an 
orbit  modification  will  be  accomplished  are  fuel  and  time. 
If  it   is  desired   to  achieve   the  new   orbit  as   rapidly  as 
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possitle,  then  a  large  amount  of  fuel  will  be  required  to 
provide  the  necessarily  large  force  and  acceleration-  Such 
a  time  constraint  may  exist  in  the  case  of  rescue  missions 
or  anti-satellite  weapons.  If  time  is  not  the  deciding 
criterion  then  slower,  more  fuel  efficient  transfers  are 
possitle.  This  has  the  major  advantage  of  allowing  more  of 
the  satellite's  weight  to  be  in  instrumentation  since  less 
weight  need  be  allocated  to  fuel,  A  fuel  efficient  transfer 
will  be  used  for  this  thesis  to  maximize  operational 
payload. 

Figure  2.4  depicts  a  satellite  at  position  one  in  a 
counter-clockwise  circular      orbit   about      the   Earth-  It   is 

desired  to  discharge  a  payload  from  the  satellite  so  that 
upon  completion  of  the  maneuver,  the  payload  will  be  in 
position  two  and  the  satellite  will  return  to  position  one- 
The  positions  are  separated  by  the  angle  u  measured  in 
radians.  This  same  type  Df  maneuver  could  be  used  if  posi- 
tion one  were  to  maneuver  in  crder  to  rendezvous  with  a 
second   satellite      already   at      position    two.  To   accomplish 

either  task,  the  satellite  at  position  one  will  accelerate 
into      the      elliptical      orDit      shown-  While      position      two 

continues  on  its  circular  orbit  and  may  pass  the  original 
position   one    several    times,  the   accelerated    satellite    will 

travel  the  elliptical  path  and  rendezvous  will  occur  when 
both  satellite  one  and  two  meet  at  position  one-  Note  that 
tne  satellite  at  position  one  can  not  merely  "speed  up"  to 
catch  position  two  because  by  Equation  2-7  they  would  be  at 
different  distances  from  the  Earth.  It  is  necessary  to  know 
how  much  position  one's  velocity  must  be  increased  to  place 
it  into  the  correct  elliptical  orbit  so  rendezvous  will 
occur . 

Since  the  satellites  are  to  meet  at  position  one, 
the  time  for  one  elliptical  orbit  must  be  equal  to  the  time 
for  an  integer  number  of  circular  orbits  plus  a  portion  of  a 
circular        orbit      to        account      for         the      initial        angular 
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Figure  2-4    Circular  Orbit  Modification 
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difference.  The  integer      number      of      orbits     (K)         can      be 

selected  to  again  either  minimize  time  or  fuel  expended  with 
more    circular   orbits    indicating   less   fuel  required-      That   is 

KTc     +    Tc(1-u/27r)    =    Te  (2.13) 

Substituting  Equations  2-10  and  2.11  into  Equation  2.13 
gives 

K27rRc^/2     ^     27rRc  ^/2  (1-u/2  7r)      =     Iwd^/^  (2.14) 

SQfiT(GM)  SQRT(GM)  SQRT  (GM) 

The  circular  orbit  radius  is  now  being  represented  by  R^ 
rather  than  just  r  to  assist  in  equation  development.  The 
common   term    of   27r/S2RT(GM)    can    be   cancelled    which    leaves 

KRc   3/2     +     £^3/2   (  1-u/27r)       =     d^  /  2  (2.15) 

Rearranging  and  solving  for  a  gives 

a     =    R^  [K  +  1-U/2  7r]2/3  (^,  1o) 

Equation    2.16    gives    an    expression    for    a    which   is   tae 
semi-major      axis      of       the   required      elliptical      orbit-  An 

expression  is  needed  that  relates  the  eccentricity  of  the 
transfer  ellipse  to  K.  This  can  be  done  by  using  Equation 
2.3  and  realizing  that  at  perigee  H^  =  fip.  ;iax.ing  the 
substitution    yields 

a    =    a  (1-e)[  K+1-U/2  7r]2/3  (2.17) 

Cancelling    the   a    and    rearranging    gives 

1     =    [  K  +  1-U/27r]2/3     -     e[  K+1-U/27r]2/3  (2.18) 

Solving    for   e    gives    Equation    2.19 
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e    =   [  K  +  1-u/27r]2/3    -    1  (2.19) 

[  K  +  l-U/2  7r]2/3 

Equation  2.19  now  relates  the  necessary  eccentricity 
of  tiie  transfer  ellipse  to  both  the  number  of  circular 
orbits  of  position  two  and  to  their  initial  angular  separa- 
tion. Equation  2.7  gives  the  circular  orbit  velocity  and 
Equation  2.12  describes  the  elliptical  orbit  velocity-  The 
velocity  of  the  elliptical  orbit  at  Rp  must  be  known  to  find 
the      difference      between    the      two      velocities.  Evaluating 

Equation    2. 12   at    r  =    Rp   gives 

Vp     =    SQET[GM(2/Rp- 1/a)   ]  (2.20) 

Vp    is    the      velocity    at  perigee    for      satellite  one.         Solving 
Equation   2.3    for    a   gives 

a    =    Rp  (2.21) 

1-e 

Substituting    Equation      2.21    into      Equation    2.20      and   simpli- 
fying   yields 

Vp     =    SQRT[  (GiVRp)  (1+e)  ]  (2.22) 

The  required  change  in  velocity  can  now  be  computed 
by  subtracting  V,.  from  Vp. 

Vp  -  V^  =  SQRTC  (GM/Rp)  (Ue)  ]  -  SQET[GH/Ec  1  (2.23) 

At  perigee   R^  =  Rp;   Vp,.   will  represent  the   difference  in 
vel ocity. 

Vp^  =  SQf<T(GM/E<.)  [  SQRT  (1+e)  -  1]  (2.24) 

The      square      root         of       (1    ♦   e)         can      be        approximated      by 
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1    «•- 5e   -••...    .         Also    note      that    SQRT(GM/Rc)     =   ^c  •         Thus 
Equation   2.2^    cam   be    written   as 

Vpc     =    VcCI    *    ,5e  -     1]  (2.25) 


V^e  (2,26) 


=     V^  [   (K+1-U/2  7r)  2/3    -     1]  (2.27) 

2[  K  +  1-U/27r]2/3 

Equation  2.27  is  a  key  formula  and  can  be  used  to 
determine  how  much  additional  velocity  is  required  to  ensure 
a  proper  rendezvous  transfer  orbit.  It  has  the  flexibility 
to  vary  the  number  of  circular  orbits  which  is  equivalent  to 
beinq  able  to  predetermine  the  rendezvous  time. 

Tne  use  of  this  equation  can  be  extended  to 
deployinq  numerous,  equally  spaced  satellites  from  a  single 
large  parent  vehicle.  The  parent  vehicle  would  be  in  the 
elliptical  orbit  carrying  the  satellites  to  be  deployed- 
Each  time  the  "rendezvous  point"  is  achieved  a  satellite 
would  be  ejected.  Two  satellites  would  not  be  placed  in  the 
same  position  as  the  new  rendezvous  point  is  always  u 
radians  ahead  of  the  current  one.  Because  Equation  2.27  has 
u  as  a  variaDle  parameter,  any  angular  positioning  is 
possible  and  numerous  satellites  may  be  deployed  limited 
only  by  the  total  weight  tnat  can  be  launched  from  Earth. 
3 •   Orbital  Transfer  fuel  Requirements 

The  next  parameter  that  must  be  calculated  is  how 
mucn  fuel  is  required  to  inject  the  parent  satellite  from 
the  circular  orbit  into  the  elliptical  orbit-  Comparison  of 
Equations  2.7  and  2.22  show  that  Vp  has  a  lar qer  maqnitude 
tnan  V^.  even  thouqh  both  are  at  a  distance  Rp  from  the 
center  of  the  Earth,    Thus,   a   transfer  from  a  circular  to 
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elliptical  orbit  requires  the  satellite  to  accelerate  in  the 
direction  of  travel.  Similarly,  when  the  parent  satellite 
ejects  a  payloai,  that  payload  must  burn  fuel  and  slow  down 
to  Vc  or  it  would  continue  on  the  parent's  elliptical  path- 
The  amount  of  fuel  required  for  either  maneuver  is  critical 
since  not  only  the  fuel  but  also  its  container  subtract  from 
the    weight   that  could   be   used   for  instrumentation. 

Rocket  fuel  energy  is  classified  by  its  specific 
impulse  (Isp)  and  is  a  measure  of  the  thrust  developed  as  a 
function  of  fuel  mass  flow  rate.  The  mass  flow  rate  can  be 
in  terms  of  gallons  per  minute  for  a  liquid  fuel  or  pounds 
per  second  for  a  burning  solid  fuel.  The  accepted  English 
units   for   Ijp    are    expressed    as    seconds.       This  comes    from 

Ijn       =      pounds    of    thrust  =      seconds  (2.28) 


sp 


pounds   of    fuel/second 


Typical  values  of  Ijp  range  from  275-325  seconds.  Knowing 
the  fuel  Isp/  the  force  of  acceleration  or  thrust  (Fj)  can 
be    found   by 

F^     =    Isp(dW/dt)  (2.29) 

W  is  the  propellent  weight  flow  rate. 

The  rocket  thrust  is  generated  from  burning  fuel 
that  exits  the  rocket  at  a  velocity  V^  with  units  of  feet 
pec  second.  The  fuel  Ijp  and  exhaust  velocity  can  oe 
related  by 

Vx  =  gisp  (2.30) 

where  g  is  the  acceleration  constant  equal  to  32-2  ft/sec^. 
The  thrust  can  also  te  related  to  the  exhaust  velocity  and 
fuel    flow    rate   by 

F^     =    -(dm/dt)  Vx  (2,31) 
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The    minus   sign  is    necessary   because  the   exhaust  is   traveling 

in    one   direction   and    the   resultant      force  is   in  the    opposite 

direction.  By  Newton*s  laws,  the  rocket  must  react  to  the 
fuel    thrust   by 

F    =    Ms^  (2.32) 

ilj  is  the  mass  of  the  rocket  being  accelerated  and  a  is  the 
resultant  acceleration.  The  acceleration  can  be  written  in 
terms   of   the   rocket    velocity    (V^ )    as 

a    =    dV^/dt  (2-33) 

Substituting    Eguation   2.33    into   Equation    2.32   gives 

F    =    Ms  (dVr /dt)  (2-34) 

Fj     must    equal      F-  Setting    Equation    2.34      equal    to    Equation 

2.  3  1    gives 

Ms  (dVr  /dt)     =   -(dm/dt)Vx  (2-35) 

Simplifying  and    rearranging    terras  gives 

dVr     =    -Vx  (dm/Ms)  (2.36) 

The  left  side  of  this  eguation  can  now  be  integrated  from 
tne       initial      rocket      velocity     (Vj  )  to      the      final      rocket 

velocity  (Vj ) .  The  right  side  can  oe  integrated  from  the 
initial   rocket   weight    (Wj  )     to    the   final   rocket    weight    (>Jj). 

/dV,       =      -V^    fdMs  (2.37) 

Equation    2.38    shows    the   evaluation    result, 
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Vf   -Vj        =       -Vxlli(Wf  /Hj) 


(2.  38) 


Dividing  by  -V^    and    using  both    sides   as   a   power    of    the   expo 
nential   e  gives 


expt-CVj-Vj  )/Vx  ]      =      Wj/Wj 


{2.39) 


W,       =      Wi     exp[-(Vj-V;)/Vx  ] 


(2.40) 


The  initial  weight  must  be  the  sum  of  the  final  weight  plus 
the  weight  of  any  fuel  or  propellent  that  was  burned  (Wp)  to 
produce    the  change  in    velocity,    or 


Wj      =    Hj  -Wp 


(2.41) 


Substituting    Eguation   2.41    into   Equation    2.40   gives 


Rj      Wp 


Wj     exp[-(Vf-Vj  )/Vx] 


(2.42) 


Solving    for    Wp   gives 


ffp      =      W,  [1    -    exp-(Vj-Vi  )/Vx  ] 


(2.43) 


Replacing    V^    with    gljp    from    Equation    2.30    gives 


Wp       =       Wjfl    -    exp-(Vj-Vi  )/(gIsp)  ] 


(2.44) 


The    quantity    V^-Vj        is  the    same   change      in    velocity    required 
by    Equation    2.27    for    orbital    transfer. 


Vpc     =    Vj-Vi 


(2.45) 


Equation    2.44    then   becomes 
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Wp    =   WiCI   -   exp-(Vpc /gisp)  ]  (2.46) 

Equation  2-27  can  be  used  to  determine  the  amount  of 
velocity  change  that  is  needed  for  an  orbital  transfer  and 
Equation  2.46  is  used  to  determine  the  percentage  of  initial 
satellite  weight  that  must  be  reserved  for  the  fuel  to 
accomplish   the  same   transfer. 
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III.    SATELLITE   DESIGN 

A.  GENEBAL    DESIGN 

The  satellite  will  be  placed  in  a  circular  low  Earth 
orbit  (LEO)  of  322  km.  An  orbit  close  to  the  Earth  will  be 
used  because  of  the  anticipated  low  power  level  of  the 
earthquake  electromagnetic  signal.  The  design  of  the  satel- 
lite will  be  based  on  using  the  space  shuttle  for  achieving 
the  initial  circular  orbit.  Because  of  the  large  cargo  bay 
on  the  shuttle  (60  ft  X  15  ft) ,  a  group  of  four  satellites 
arranged  in  a  circle  and  mounted  on  a  central  support  bus 
will     comprise   the      launch      package.  When  deployed,         the 

satellites    will    be   equally    spaced   by    90    degrees. 

Primary  electrical  power  will  be  provided  by  solar  cells 
with  tatteries  supplying  power  during  the  time  the  satellite 
is  in  the  Earth's  shadow.  Spin  stabilization  of  the  satel- 
lite body  With  hydrazine  thrusters  for  station  keeping  will 
be  used  to  maintain  a  de-spun  antenna  pointing  at  the  Eartn 
with  attitude  sensing  to  be  done  by  Earth  limb  sensors. 
Each  satellite  will  weigh  approximately  500  pounds  and  have 
a   cylindrically   shaped   main    body. 

B.  AMTENHA    DESIGH 

The  satellite  will  be  designed  to  detect  electromagnetic 
energy  in  the  30  MHz  region.  This  frequency  is  used  because 
it  is  designated  for  astronomical  observations  and  manmade 
interference  should  be  at  a  minimum.  Additionally  30  flHz  is 
well  above  the  cutoff  frequency  (Fp)  for  propagation  through 
the  ionosphere  given  by 

Fp  =  9SQRI(Ne)  (3,1) 

where  N©   is  the   ionospheric  plasma   electron  density   with 
units   of  electrons/meter^.     a   nominal   value  for   Fp   is 
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10    MHz.  Any    frequency      higher   than      Fp    will      successfully 

penetrate  into  outer  space  with  lower  frequencies  being 
reflected   back   to   Earth, 

A  Yagi  antenna  will  be  used  as  the  main  sensor  because 
of  the  size  and  space  constraints  imposed  by  a  satellite 
design.  fieference  8  will  be  used  as  a  basis  for  the  design 
of  the  antenna.  Table  I  [Ref.  8:  p.  7]  is  used  once  the 
initial  design  criteria  is  established.  This  initial  data 
consists  of  the  frequency  ot  operation  (30  MHz) ,  antenna 
gain  desired  (9.2  dB ) ,  d/A  ratio  where  d  is  the  diameter  of 
the  parasitic  elements  (2  cm),  and  the  D/a  ratio  where  D  is 
the  diameter  of  the  support  boom  (U  cm)  .  The  wavelength  A 
is    given    by 

A  =  c  (3.  2) 

f 

The         speed      of        light       (c)  is      a         constant      equal        to 

3  X  108  meters/second  and  f  is  the  operating  freguency- 
Thus  a  frequency  of  30  ;iHz  gives  a  wavelength  of  10  meters- 
Following  the  steps  in  Reference  8,  Figure  3.1  [Sef.  8:p. 
20]  gives  the  length  of  the  first  and  third  director  as 
4.U4  meters,  4.5  meters  for  the  length  of  the  second 
director  and  4.9  metecs  for  the  length  of  the  reflector. 
The  element  spacing  is  2  meters.  Deployment  of  the  antenna 
in  outer  space  can  be  achieved  by  using  a  bistem  whicn  is  a 
deployable  boom  rolled  up  on  a  drum  in  its  flattened  config- 
uration. Because  of  its  small  size  when  stowed,  total 
storage  space  required  will  be  approximately  one-half  a 
cubic  foot.  Since  the  antenna  will  be  operating  in  a 
weightless  environment,  strong  structural  members  are  not 
required.  The  total  weight  of  the  antenna  would  not  exceed 
30  pounds  since  the  material  of  construction  could  be  a 
light-weight  aluminum  based  alloy.  The  gain  of  9.2  dB  is 
referenced      to      a      dipoie   antenna.  For      reference      to      an 

isotropic   radiator,    2.16    dB    of    gain    must    be    added. 
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TABLE   I 
YAGI    ELEMENT    LENGTHS 


LENGTH   Of    tACI    IN  WAVELENGTHS 
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I  Ith 
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0.386 
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0.20 

0.20 

0.25 

0.20 

o.;o 

0.308 
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C.       FUEL    EEQDIBEHEHTS 
1.      Satellite  Fuel 

The  on-station  weight  of  the  satellite  is  500  pounds 
and  this  weight  must  be  decelerated  from  the  parent' s  ellip- 
tical velocity.  Equation  2.27  will  be  used  to  determine  the 
necessary  velocity  change  and  Equation  2.7  will  be  used  to 
find        Vj  •  GM        evaluates        as        a      constant        equal        to 

3.99  X  101*  m3/sec2  and  r  is  the  sum  of  the  radius  of  the 
Earth  {6370  Km)  and  the  satellite  altitude  above  the 
Earth    (322    Km)  . 

V^      =    SQRT(3.99    X    10i*/6692    X    103)  (3.3) 


=    7,722    m/sec  (3,4) 


=   25^335    ft/sec  (3.5) 

Equation   2.27    can    now    be   evaluated    if  K    is    set    equal    to   one 
and    u    is   set   equal   to     -a /2    radians. 

Vp,     =    25,335[  (2-. 25)2/3    -     i]  (3.6) 
2[  (2-. 25)  2/3     -    1  ] 

=  3,945    ft/sec  (3.7) 

Equation   2.46    can    be    used   to   find    the  fuel   fraction    given    by 

Wp/Wj     =    1    -    exp-[  3945/(32)  (300)  ]  (3.8) 

=    .337  (3.9) 
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Figure    3.1         Yagi    Antenna    Design   Curves 
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The  initial  weight  is  composed  of  the  en- orbit 
satellite  weight  of  500  pounds  (Hs ) #  ^p  and  the  weight  of 
the    case    (W^)     that  contains    the  fuel.      That    is 

Wi      =    iis+Wp  +  Wc  (3.  10) 

For  an  efficiently  designed  rocket  motor,  the  majority  of 
the  weight  is  due  to  fuel  with  the  remaining  weight  due  to 
the    case.      This  ratio    (f)    is   given    by 

f    =_Wp (3,11) 

Wp*Wc 

A  typical  value  of  f  is  0.93.  Using  this  definition  of  f, 
Eguation   3.10    can   be    rewritten   as 

Wj      =    Hs+Wp/f  (3-  12) 

Substituting  Hp  =  .337Wi  from  Equation  3,9  and  0-93  for  f 
into    Equation   3.  12   gives 

Hj      =    Ws+.337Wi  /.93  (3.  13) 

Rearranging      terms        to      solve        for      Wj  and      substituting 

500    pounds   for   Hs     gives 

Hj     =  _500 (3.  U) 

1    -    .337/, 93 

=    784    lb  (3. 15) 

This  can  be  used  to  find  Wp  from  Equation  3,9 

Wp  =  (784)  (.337)  (3,  16) 
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Wp     =    264    lb  (3. 17) 

The      ueight      oi  the      case        can      be     found      by      rearranging 
Equation   3.11 

He      =    Wp(l-f)  a.  IB) 


=    264  (1--93)  (3.  19) 

.93 

=    20   lb  (3.20) 

2  .      Parent  Satellite    Fuel 

The  amount  of  fuel  required  by  the  parent  to  inject 
itself  into  the  elliptical  orbit  can  now  be  calculated  in  a 
similar  manner-  The  central  support  bus  weight  will  be  set 
at    500    pounds-  This  relatively      large    weight      is    required 

because  the  bus  must  be  strong  enough  to  rigidly  hold  the 
four    satellites      during   the    sh uttle- launch.  Thus,      enough 

fuel  must  be  allowed  to  accelerate  the  bus  and  three  satel- 
lites into  an  elliptical  orbit  since  one  satellite  is 
deployed  prior  to  the  parent's  acceleration.  The  net  weight 
to  be  boosted  is  then  2,852  pounds.  Similar  to  Equation 
3.14,    the   initial    total    weight    is 

Wj      =   2852 (3-21) 

1    -    .337/. 93 

=    4,473    lb  (3.  22) 

The  fuel  weight  and  case  weight  for  the  parent  can  also  be 
found    in   the    same    manner    as    for   the   individual    satellite. 

Wp    =     (4473)  (.337)  (3.  23) 
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Hp     =    1,507    lb  (3-24) 


Wc     =    1507(1--93)  (3.25) 

.93 


=    114    lb  (3-26) 


D.       ELECTRICAL   POHEfi 

The  primary  source  of  electrical  power  will  be  provided 
by  converting  light  energy  into  electrical  energy  by  use  of 
solar  cells-  Because  the  available  power  from  an  individual 
cell  is  small  (on  the  order  of  0.06  watts),  a  large  number 
of  cells  must  be  used  that  collects  solar  energy  over  a 
large  incident  area.  In  order  to  provide  the  desired  total 
power  at  useful  voltages,  solar  cell  configurations  are 
normally  used   in    a   series-parallel    connection. 

The  satellite  will  not  be  in  sunlight  during  the 
Earth-Sun  eclipse  periods  of  the  orbit-  During  these  times, 
rechargeable  batteries  will  provide  the  electrical  power. 
Additionaly,  the  batteries  can  provide  back-up  power  during 
transients   that   exceed  the   solar    panel    capacity. 

In  order  to  provide  a  useful  power  system,  these  two 
power  sources  must  be  integrated  with  a  power  control  and 
power      conversion   system.  Functions    that      must    be      accom- 

plished include  battery  charging  during  sunlight  periods, 
bus  voltage  monitoring  and  control  and  converting  battery 
power    to  electrical    power   when    required. 

Figure  3.2  [  Ref.  9:p.  1-7.5]  shows  a  block  diagram  of  a 
satellite  power  system.  While  only  one  of  each  component  is 
shown,  normal  installations  consist  of  multiple  components 
for  reliability  reasons.  The  power  from  the  main  solar 
panel  supplies  the  spacecraft  systems  through  the  primary 
power   distribution   bus.      The   output    voltage   of    the    panel    can 
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Figure  3.2   Satellite  Power  System 
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vary  as  a  function  of  load,  temperature  and  other  factors. 
Bus  voltage  limiters  prevent  the  voltage  swings  from 
exceeding  spacecraft  design  ranges  by  acting  as  variable 
shunt  loads  during  a  rising  voltage  condition.  The  battery 
discharge  regulator  turns  on  when  a  pre-determined  bus  low 
voltage   is      sensed  and      provides  a      regulated   output.  The 

battery  terminal  voltage  is  higher  than  the  required  bus 
voltage  and  in  order  to  charge  the  battery  a  still  higher 
level  of  voltage  is  needed.  To  provide  this  higher  voltage 
level,  the  additional  solar  panel  is  series  connected  to  tne 
main  panel.  The  size  of  the  battery  charging  panel  is 
designed  so  that  the  charge  control  functions  provided  in 
the  controller  need  only  be  an  on-off  and  rate  change 
switch.  The  voltage  and  current  sensors  provide  means  to 
monitor  the  operational  state  of  the  power  system  and  for 
load  managing  during  the  orbit. 
1  -      Solar    Pa n el    Design 

As  shown  in  Figure  3.3  [Ref.  9:p.  1-7.11]  individual 
solar  cells  are  arranged  in  a  series-parallel  array  to 
provide      the   required      power.  The      three    strings      are      in 

parallel  to  minimize  power  loss  should  a  single  cell  fail. 
These  arrays  are  mounted  on  the  satellite  with  the  series 
string  length  parallel  to  the  satellite  spin  axis.  All  the 
arrays  are  then  connected  in  parallel  to  produce  the  desired 
power.  The  surface  of  the  solar  cell  is  considered  to  oe  a 
Lambertian  surface  and  the  power  produced  by  the  cell 
decreases  by  a  cosine  function  as  the  sun's  incident  angle 
changes  from  normal  as  shown  in  Figure  3.4  [Bef-  9:p. 
1-7.  11]. 

Figure  3.5  [Bef-  9:p.  2-4.17]  shows  how  the  sun's 
incident  angle  varies  during  one  year.  When  the  solar  cells 
are  in  the  sunlight,  only  those  positioned  in  the  center  are 
illuminated  by  normal  rays.  Going  left  or  right  from  center 
causes  the  ceil  output  to  drop  by  the  cosine  function.  The 
net    result    is    an    output   from   the   half -cylinder    that    is    equal 
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Figure   3-3        Solar  Cell   Array 

to  the  output  from  a  flat  plate  whose  dimensions  are  the 
diameter  and  height  of  the  cylinder.  Additionally,  the  sun 
line  varying  23-5  degrees  as  shown  in  Figure  3-5  causes  a 
further   decrease    in    solar   panel   output   power. 

Solar  cell  performance  also  degrades  with  time  from 
exposure  to  electron  and  proton  fields  trapped  in  the 
Earth's    magnetic      field  and      solar    flares.  The  extent      of 
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SUN 


Figure  S.'l        Cylindrical   Panel   Installation 

degradation      is   shown      is      Figure      3,6    [  Ref .    9:  p.  1-7.13]. 

Damage  is  a  result  of  charged  particles  creating  defects  in 
the  n  and  p  layers  of  the  cell.  The  net  result  is  shown  in 
Figure  3.7.  For  any  set  of  four  sequential  peaks,  the  maxi- 
fliums  correspond  to  the  autumnal  and  vernal  equinox  and  the 
minimuffls  correspond  to  summer  and  winter  solstice  when  sun 
incidence   is   offset    by   23-5    degrees. 
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Figure    3-5        Solar   Incident    Angle 

The  satellite  that  is  being  proposed  will  have  a 
diaaeter  of  5-7  feet  and  a  height  of  10  feet.  This  gives  a 
projected  area  of  57  square  feet.  At  beginning-of-lif e 
(BOL)  solar  cells  produce  approximately  11  watts  of  power 
per    square    foot    of   area.      Thus    the    57    square   feet    will    yield 
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Figure    3.6        Solar  Cell  Radiation   Damage 

627  watts  of  power.  This  initial  power  level  will  decrease 
because   of      the   factors   mentioned    previously.  From  Figure 

3.6,  after  seven  years  on-station,  the  power  level  will  have 
fallen  to  approximately  72  percent  of  the  BOL  value-  This 
leaves  451  watts  of  available  power.  A  second  correction 
factor  can  be  obtained  from  Figure  3.7.  Daring  any  given 
year  the  lowest  power  output  is  during  summer  solstice.  The 
highest  power  output  is  the  immediately  preceding  vernal 
equinox-  Using  the  time  during  the  third  to  fourth  years  on 
the    graph,    the  correction  factor  is 


553/625    =    .8848 


(3.27) 


The       power    that      will      be      available   at      end-of-life      during 
summer    solstice   will    then      be    approximately    400    watts.         The 
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electrical  power  budget  must  be  based  on  this  EOL  value- 
The  power-to-mass  ratio  at  EOL  for  solar  ceils  is  approxi- 
mately 20  watts/kg  for  normal  incidence  [ Hef .  10:p.  73]. 
This  ratio  must  be  decreased  because  the  satellite  does  not 
have  flat  plate  collectors.  The  amount  of  decrease  will  be 
the  ratio  of  the  area  of  a  flat  plate  collector  to  the  area 
of  the  illuminated  half  cylinder.  The  area  of  an  equivalent 
flat    collector   is   57    square    feet.  The   surface   area   of   the 

half    cylinder    (A^  )    is 

Ac    =    (-5)  (27rrh)  (3.28) 


=    (tt)  (5.7/2)  (10)  (3.29) 


=    89.5    ft2  (3-30) 

The    power-to-mass   ratio   of    the    cylinder       (P^  )    is   then 

Pc     =    (20)  (57/89.5)  (3.31) 


=    12.73   watts/Jcg  (3.  32) 

The    mass   of    the   solar   cells     (M^)    is    then 

Mc     =    200/12.73  (3.33) 

=    15.7    kg  (3.  34) 

This    is      the    mass      of   only      those   cells  that   are      receiving 
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sunlight.  The  total  mass  is  twice  as  much.  Converting  kg 
to  pounds  will  give  a  total  weight  of  the  solar  cells  as 
69    pounds- 

2-      Batter  y  Requirements 

Batteries  that  are  used  in  a  weightless  environment 
can  not  be  of  the  "flooded"  electrolyte  type  such  as  lead- 
acid  cells.  fiather,  the  electrolyte  must  be  absorbed  in  the 
separator   material      to  assure      even   distribution.  The   two 

primary  types  of  batteries  that  meet  this  condition  are  the 
sealed  silver-zinc  and   the    nickel-cadmium. 

Table  II  [Ref.  9:p.  1-7-17]  shows  the  major  charac- 
teristics of  concern  for  both  battery  types.  Nickel-cadmium 
batteries  will  De  used  because  of  their  longer  life  expec- 
tancy and  tolerance  to  overcharge  allowing  a  simpler  charge 
control  circuit.  The  major  disadvantage  will  be  an  increase 
in    weight   due    to    the    lower   energy  density. 

The  satellite  will  be  spending  50  percent  of  its 
time  in  the  Earth's  shadow  during  which  tne  battery  must 
supply  the  electrical  power.  During  the  time  sunlight  is 
available,  the  battery  must  be  recharged.  Since  equal 
amounts  of  time  are  spent  in  light  and  shade,  half  of  the 
solar  panels  must  be  dedicated  to  battery  charging  and  the 
remaining  half  to  satellite  bus  power.  The  time  the  satel- 
lite is  in  aarkness  (Tj)  can  be  found  from  Equation  2.10  to 
be 

Tj     =    .5Tc  (3.35) 


=     (tt)   (6692    X     103)  3/2/SQRT  (3.99    X    IQi*)  (3.36) 


Tj     =    45,4    min  (3.37) 
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TABLE   II 


BATIEBY    CHARACTERISTICS 
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The  satellite  system  will  be  drawing  200  watts  of  power 
during  this  time  and  the  watt-hours  required  from  the 
battery    (Pj^)    will   be 

Pjj    =    (200)  (45,4/60)  (3.38) 


=    151.3    watt-hours  (3.39) 

Figure      3.8    [Eef-    9 :  p.  1-7-21]    shows      tnat    as      the 

depth  of  discharge  (DOD)  decreases  the  lifetime  of  the 
battery  increases.  In  order  to  have  a  long  lifetime,  a  high 
capacity  battery  system  is  desirable  but  the  penalty  is 
additional  weight.  From  Reference  10,  the  battery  system  on 
the  Intelsat  V  can  supply  34  amp-hours  at  27  volts  and  has  a 
mass  of  32  kg  [Ref-  10:p.  73]-  This  gives  a  battery  energy 
density  of  13  watts/pound  which  is  slighty  better  than  shown 
in  Table  II.  To  provide  a  low  DOD  but  not  require  excessive 
weight,  the  battery  system  for  the  proposed  satellite  will 
be    capable   of    supplying    1200    watt-hours.  This  gives    a    DOD 

of 

DOD    =     15  1-3/1200  (3.40) 


=    12. o    percent  (3-41) 

From  Figure  3.3  this  gives  approximately  23,000  cycles  for 
the  battery.  The  time  for  one  cycle  is  90.8  minutes.  There 
are  5,789  cycles  completed  in  one  year.  Thus  the  battery 
has  a  lifetime  of  approximately  four  years  and  will  weigh 
about    92    pounds. 
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Figure    3,8 


Nickel-CadBium    Battery  Cycle  Life 
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E.       SYSTEM    WEIGHT   AND    POSEfi    BODGET 

The  satellite  has  an  on-orbit  weight  of  500  pounds. 
Subtracting  30  pounds  for  the  antenna,  69  pounds  for  the 
solar  cells  and  the  battery  weight  of  92  pounds  leaves 
309  pounds  for  the  remaioing  electronics,  structure  and 
control  systems.  Using  the  weight  and  power  ratios  in 
Reference    10   gives  the  results   shown  in    Table   III. 

I 

I 

TABLE    III 
WEIGHT    AND    POWER    BUDGET 


System 

Weight  (lb) 

Power (watts) 

Antenna 

30 



Solar 

Cells 

69 



Battery 

92 



Telemetry  , 
Command 

16 

7.9 

Attitude 
Control 

43 

9.6 
1.0 

Propulsion 

24 

Structure 

86 



Electrical 
Integration 

25 

20.3 

Signal 
Process in  g 

1  15 

161.2 

... 

Of  the  161-2  watts  available  for  signal  processing, 
approximately  30  percent  would  be  available  for  use  to 
transmit  data  back  to  central  Earth  collection  stations. 
The  signal  strength  at  the  surface  of  the  Earth  can  be 
calculated    by 

P,      =    ?^  (G^GJ  A2  (3.42) 

(47r)2r2 
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Pi  is  the  transmitted  power,  G^  and  G^  are  the  gains  of  the 
transmitting  and  receiving  antennas  and  r  is  the  distance 
between  the  satellite  and  Earth.  Jsing  the  designed  antenna 
as  the  transmitting  antenna,  the  gain  is  11.36  dB  with 
reference  to  an  isotropic  radiator-  An  Earth  based  antenna 
does  not  have  the  same  weight  constraints  and  could  be 
designed  for  a  gain  of  150-  Using  these  values  ,  tne 
received  signal   power  would    be    0-63    microwatts, 

F.       STABILITY 

The  satellite  is  designed  to  spin  about  its  longitudinal 
axis.  If  it  were  to  tumble  about  its  transverse  diameter, 
it  would  be  totally  inoperative-  An  analysis  of  the  momeat 
of  inertia  about  both  axes  is  required  to  ensure  proper 
placement  of  the  internal  components.  From  Reference  11  the 
moment  of  inertia  for  a  hollow  circular  cylinder  of  length  h 
about  the  longitudinal  axis  is 

Ig  =  {m/2)  (EjZ+R^S)  (3.U3) 

Rj  and  £q  are  the  inner  and  outer  radii  and  m  is  the  mass  oz 
the  cylinder.  The  moment  of  inertia  about  the  transverse 
diameter  is 

I^     =   m[Ei2  +  R^,2    +    h2]  (3.  44) 

4  12 

for  stability,  Ig  must  be  greater  than  I^,  The  longitudinal 
moment  of  inertia  due  to  the  solar  cells  with  Ej=Ro  yields 
560-4   ft^-lb-  The      corresponding         transverse      moment      of 

inertia  is  655,2  ft^-lb.  In  order  to  increase  Ig  more  than 
Ij,  weight  must  be  distributed  in  a  dj sK-like  manner  at  the 
center  of  the  satellite.  The  weight  of  the  antenna,  solar 
cells  and  structure  can  not  be  centralized.  The  remaining 
weight  is  315  pounds-  This  weight  with  Rj  equal  to  zero 
adds    1279.3    ft2-lb   in    the   longitudinal    direction   for    a    total 
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of    1839.7  ft2-lb.  It  is      dow   necessary      to   determine      the 

maximum   height    (d)    of    the   disk.      The   total   transverse   moment 
of    inertia   must  be  less  than    1839.7    ft2-ib.      That  is 

1839.7   >     (315/2)[  .5  (5.7/2)  2+d2/6]+855.2  (3.45) 


6.25-4.06    >    d2/6  <3.  46) 


d    <    3.6   feet  (3.  47) 
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17-     DESCRIPTION    OF    BESEABCH 

A.       DESIGN    OVEfiVIBM 

The  earthquake  monitoring  system  is  installed  three 
miles  east  of  San  Juan  Batista,  Ca.,  and  seven  miles  south 
of  Hollister.  It  is  located  in  the  Hollister  Hills  State 
Vehicular  Recreation  Area  (Upper  Ranch)  south  of  the  MX 
track       (see   Figures      4-1      and    4.2)-  The      site   was      chosen 

because  it  is  both  historically  and  recently  an  active 
earthquake  area  and  is  far  enough  from  industrial  areas  to 
minimize  manmade  noise  signals  yet  near  possible  earthquake 
sources  to  generate  electromagnetic  energy.  Exact  location 
is  1000  feet  from  the  San  Andreas  fault  at  121  o  23.5'  W  and 
360    45.5'    N. 

Two  frequencies  (38.45  MHz  and  150.75  MHz)  are  being 
monitored-  These  frequencies,  each  with  a  50  KHz  bandwidth, 
have  been  assigned  to  this  project  by  the  Army  Frequency 
Coordinator  for  the  Western  United  States  for  a  renewable 
one      year    period-  The      antenna      complex   consists      of      two 

20  foot  towers  and  mounted  antennas  as  shown  in  Appendix  D- 
Each  tower  has  five  mounted  antennas.  Eight  Yagi  antennas 
sense  horizontal  and  vertical  polarization  at  each  frequency 
both    up    and    down    the    fault.  Additionally,      an    omni    antenna 

at    each    frequency    is    vertically   mounted. 

Eacii  antenna  signal  is  individually  cabled  to  a  modular 
AM  receiver.  Output  from  the  receiver  is  sent  to  an  inter- 
face amplifier  where  it  is  joined  with  a  system  timing 
signal-  Tne   composite      signal      is    then      sent      to    a      cnart 

recorder  and  permanently  traced  onto  pressure  sensitive 
paper.  Figure  4.3  shows  a  block  diagram  for  each  channel 
With  a  single  timer  signal  being  common  to  all  ten  channels- 
All  components  are  powered  by  12  V  DC.  Collected  data  will 
be  correlated  to  earthquake  data  provided  by  the  Geological 
Survey   Office,    Menlo    Park,    Ca. 
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Figure   4-1         Site   Location 

B-        SYSTEH    COHPONENTS 
1  -      Eeceivers 

The  antenna  signal  is  coupled  into  the  RF  amplifier 
by  a  50  ohm  RCA  jack.  The  oscillator  frequency  is  deter- 
mined by  the  crystal  Y1.  For  38,45  MHz,  the  crystal 
frequency  is  49.15  MHz  and  for  150.75  MHz  the  crystal 
frequency  is  46.68333  MHz.  Figure  4.4  shows  the  schematic 
diagram  for  the  150.75  MHz  receiver.  The  38.75  MHz  receiver 
does    not   have   the   Q4    multiplier  circuit. 
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Figure   4.2         Site   Location 
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The  10.7  MHz  IF  signal  is  sent  to  the  IF-20  IF/Audio 
assemi3ly.  The  signal  is  mixed  with  the  10-245  MHz  oscil- 
lator frequency  and  down-converted  to  455  KHz-  This  signal 
is   then    amplified   and   detected-  The   "S"    meter  output    from 

an    integral  log   detector   is    used      as    the    input    to   the  inter- 
face   amplifier.         Figure    4.  5   is     a    schematic   diagram   for   the 
IF/detector  circuit-         Appendix  B   contains   measured   receiver 
sensitivity   curves- 
2  -      Timer 

In  order  to  determine  the  time  of  data  reception, 
the  timer  generates  a  voltage  signal  every  15  minutes, 
hourly  and  every  24  hours.  The  duration  of  this  signal  was 
modified  from  the  original  design  to  new  values  of 
14  seconds  for  the  15  minute  pulse,  42  seconds  for  the 
hourly  pulse  and  98  seconds  for  the  24  hour  pulse-  Accurate 
time  is  required  to  allow  comparison  with  earthquake  data. 
A  21  stage  counter  timer  (RCA  CD4045A)  is  crystal  controlled 
and  establishes  a  one  second  time  reference  for  the  entire 
system.  Synchronous  down  counters  (RCA  CD40103B)  then  use 
this  reference  to  generate  signals  at  the  appropriate  times- 
The  reset  button  is  used  to  initialize  or  reset  the  zero 
time  mark-  A  clock-running  LED  is  used  as  a  visual  indica- 
tion of  proper  operation  and  a  pen-mark  LED  lights  for  each 
occurance  ol  a  time  pulse  and  stays  lit  for  its  duration- 
The  time  pulse  shows  on  the  chart  paper  as  a  momentary 
increase  in  signal  level  and  is  above  the  local  average 
recorded      signal    strength-  Appendix    C    contains      examples. 

Figure    4.6   is    the   schematic    diagram- 
3 .      Interface    Amplifier 

The  "S"  meter  output  from  the  receiver  and  the  time 
pulse  from  the  timer  are  integrated  into  one  signal  in  the 
interface  amplifier-  This  is  accomplished  by  using  opera- 
tional amplifiers  (LM  148/348)-  On  Figure  4.7,  R3  is  used 
as  a  zero  adjust  to  establish  a  zero-signal  reference 
voltage-        ElO    adjusts  the    current    to    the   chart    recorders   so 
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Figure  4.3        Monitoring   System   Block  Diagram 
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the  meter  stylus  can  be  set  to  its  zero-signal  position.   R5 
is   then  adjusted  for  full   scale  meter   deflection  at   the 
desired  RF  signal  strength  and  acts  as  a  gain  control. 
4 .   Chart  Recorders 

The  Eustrak  recorder  stylus  is  gear  driven  and 
strikes  the  pressure  sensitive  paper  at  regular  intervals  to 
produce  an  apparently  continuous  line.  The  stylus •  position 
depends  upon  the  current  level  through  the  meter.  The  paper 
speed  is  determined  by  interchangeable  gear  trains-  Two 
styles  of  recorders  are  used  and  both  are  geared  to  allow 
one  roll  of  paper  to  last  two  weeks.  The  main  functional 
difference  between  the  styles  is  the  amount  of  current 
required  for  full  scale  deflection. 
5-   Batter  y  Backup 

Figure  4.8  shows  the  power  supply  block  diagram. 
When  60  HZ,  115  volt  power  is  available,  the  12  volt  power 
supply  is  used  to  provide  system  power.  The  relay  contacts 
are  held  open  by  the  60  HZ  power  and  the  battery  chargers 
maintain  full  charge  on  both  batteries.  The  five  diodes  in 
series  (for  both  positive  and  negative  supplies)  are 
reversed  ciased  with  no  current  flow.  In  the  event  AC  power 
is  lost,  the  relay  contacts  close  and  power  is  supplied  by 
the  batteries  through  the  three  parallel  diodes.  Three 
diodes  in  parallel  were  used  for  current  handling  and  reli- 
aoility  reasons.  When  A3  power  is  regained,  the  system 
returns  to  its  original  condition. 

It  was  important  to  maintain  power  to  the  timer 
without  interruption  or  the  time  sequence  became  scrambled. 
Tnis  is  accomplished  by  the  five  series  diodes  for  the  time 
interval  between  losing  AC  power  and  the  relay  contacts 
closing.  The  parallel  diodes  drop  the  battery  voltage  by 
approximately  0.  b  volts.  This  was  necessary  because  the 
fully  charged  batteries  produce  12.6  volts  for  several  days. 
While  not  critical  to  the  electronics  portion  of  the  system, 
the  chart  recorder  paper  drive  motors  are  DC  powered  and  the 
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Figure   4.4        150.75    MHz    7HF   Converter 
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resultant  increase  in  speed  would  result  in  running  out  of 
paper  before  the  designed  two  weeks. 
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Figure  4.5   IF-20  IF/Audio  Assembly 
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Figure  4-6        Timer    and   Reference  Generator 
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Figure  4.7    Interface  Amplifier 
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Figure    4-8        Battery  Backup   System 
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V.  BESOLTS  AND  CO NCLO SIGNS 

A.  EESOLTS 

The  system  installation  was  completed  and  began  opera- 
tion on  November  8,  1985,  The  first  set  of  strip  chart 
recordings  cover  the  time  period  of  8-22  November,  1985- 
Appendix  C  contains  samples  from  channel  seven  which  is  at 
150.75  MHz,  vertically  polarized,  and  iooicing  toward  the 
southeast.  During  the  timeframe  shown,  three  small  earth- 
quakes occured.  The  first  was  a  magnitude  2.3  at  2044  on 
November  11,  the  second  a  magnitude  2,9  at  1U11  on 
November  12  and  tne  third  a  magnitude  2,04  at  0602  on 
November  13,  No  other  channels  recorded  any  signals  during 
the  same  time  period.  The  magnitude  2,3  was  17  miles  from 
the  antenna  site  with  the  other  two  being  closer.  All  tnree 
were  in  a  southeast  direction  from  the  antennas,  Channel 
ten  (38,45  MHz,  omni)  recorded  two  groups  of  signals  witn 
the  first  occuring  during  the  three  days  preceding  the 
quakes,      Tne    second    group   y as    two    days    after    the   quakes, 

B.  CCHCIUSIONS 

While  no  firm  conclusions  can  be  reached  from  such  a 
limited  set  of  chart  recordings,  tne  electromagnetic 
activity  on  channels  seven  and  ten  lend  encouragement  to  the 
theory  of  earthquake  generated  RF  energy.  Additional 
recordings  are  required  to  confirm  this  theory  or  show  just 
a  random  coincidence  on  the  first  attempt  at  data 
collection. 

In  order  to  show  a  connection,  recordings  should  be 
continued  for  at  least  a  year.  Not  only  tne  simultaneous 
occurence  but  also  those  where  a  quake  occured  without 
recorded  RF  energy  (and  vice  versa)  would  give  valuable 
information  concerning  reliability  and  falsa  alarm  rate. 
Any  relationship  between   the  recorded  signal  level   and  the 
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magnitude  of  the  quake  or  depth  and  linear  distance  of  the 
epicenter  from  the  antennas  should  also  be  investigated. 
The  present  equipment  will  provide  an  initial  level  of 
confidence  on  which  to  base  a  future  decision  of  whether  or 
not    to  continue  the   research.  If  sufficient    evidence   does 

exist  to  justify  further  time  and  monies,  then  more  sofisti- 
cated  and  sensitive  equipment  will  be  required  to  detail  any 
further  relationships.  The  ideal  system  would  include  auto- 
mated data  collection  similar  to  those  currently  in  use  by 
the  Geological  Survey  Office.  When  sufficiently  understood, 
the  possibility  of  transferring  the  monitoring  to  a 
space-based  sensor  could  finally  be  investigated  with  the 
advantages    mentioned    earlier. 
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APPENDIX  A 
CHANNEL  CONNECTIONS 


TABLE  IV 
CHANNEL  CONNECTIONS 


Fregaency 

(Ihz) 

Channel 

V  =  Vertical 
H  =  Horizontal 
SE  =  South  East 
NW  =  North  West 


Polarization    Direction 

V  NW 

H  SE 


3 

150.75 

OMNI 

— 

4 

38.45 
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NW 

1        5 

150.75 
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SE 

1        6 

38.45 
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SE 

1        7 

150.75 
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SE 

1         8 

38.45 
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NW 
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APPENDIX    B 
RECEIVER    SENSITIVITI    CORVES 
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RECEIVER  SENSITIVITY(150.75  MHZ) 
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RECEIVER  SENSITIVITY  (38.45  MHZ) 
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APPENDIX    C 
CHART    PAPER    RECORDINGS 
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Figure  C,2   Channel  7  Recordings 
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APPENDIX    D 
PflOTOGfiAPHS   OF    DATA    COLLECTION    SISTEM 


Figure   D.I         ET1    Babka   and   Recorder  System 
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Figure   D,2        Antenna   Site    and    Protective   Building 
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Figure  D-3    Antennas  and  Towers 
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Figure  D.4    Antennas  and  Towers 
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